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ABSTRACT 



Some of the first results are reported from RISE - a new fast camera mounted on the Liverpool Telescope primarily designed to 
obtain high time resolution light curves of transiting extrasolar planets for the purpose of transit timing. A full and partial transit of 
WASP-3 are presented, and a Monte Carlo Markov Chain analysis is used to update the parameters from the discovery paper. This 



results in a planetary radius of 1.29 I .iR, 



and therefore a density of 0.82+JJ^pj, 



consistent with previous results. The inclination 



is 85.06™ [j deg, in agreement (but with a significant improvement in the precision) with the previously determined value. Central 
transit times are found to be consistent with the ephemeris given in the discovery paper. However, a new ephemeris calculated using 
the longer baseline results in r c (0) = 2454605.55915 ± 0.00023 HJD and P = 1.846835 ± 0.000002 days. 
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1. Introduction 

Transit observations of extrasolar planets are vital to our under- 
standing of planetary systems, as they allow us to determine 
the orbital inclination and planetary radius, and when coupled 
with the radial velocity (RV) method, the mass and density may 
be derived. Ground- based transit s urveys such as Su perWASP 
dPollacco et alj|2006l) and HATnet dBakos et al.ll2004l) are now 
coming of age and to date over 50 transiting planets have been 
founqj. Most of thes e systems are Hot J upiters (with the excep- 
tion of GJ-436b - see lGillon et alJ d2007)). and neither these tran- 
sit surveys nor the RV method are close to reaching the required 
precision to search for Earth-sized planets. 

However, high precision transit photometry may provide us 
with a short-cut to detecting Earth-sized planets around exist- 
ing transiting systems, as it allows us to measure the central 
transit times of these systems very accurately. This could reveal 
the presence of a third body through transit timing variations 
(TTV) by measurin g any perturbations on the transiting planets 
orbital period dMha lda-Escu dill2002t iHolman & Murrav l l2005t 
lAgol et alJ 120051: iHevl & Gladmanll2007l) . The TTV technique 
is sensitive to Earth-sized planets (or smaller) when the per- 
turber is in favourable orbital configurations - namely low-order 
mean motion resonances. The technology to take advantage of 
this technique already exists and in fact ground-based TTV ex- 
periments are the first to have probed for Earth -sized planets 
around main sequence star s (see for example, Holma rTet alJ 
2006: ISteffen & Agolll2005l) . 

The TTV method requires many light curves of high-time 
resolution coupled with precision photometry when search- 
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ing for Earth-sized exoplanets. RISE (Rapid Imager to Search 
for Exoplanets) is a new instrument designed primarily for 
this purpose. It was designed and built by scientists from 
Queens University, Belfast, Liverpool JMU and the University of 
Manchester, and is mounted on the fully robotic 2-m Liverpool 
Telescope (LT) on La Palma. It was commissioned in late 
February 2008 and is described in §2.1. 

In the months following the commissioning run some non- 
TTV targets were observed with RISE to test the instrument and 
reduction procedures further. These data have already been used 
to accurately characterise a recently disco vered plan e t from the 
SuperWASP-N camera, WASP- 14b (see Uoshi et all d2008b for 
details). In this paper we present a full transit of WASP-3 ob- 
served during this period, along with a recently obtained par- 
tial transit. W ASP-3b was also disc overed by the SuperWASP- 
N camera (see Pol lacco et al.ll2008h and was found to be a rel- 
atively massive planet closely orbiting an unevolved main se- 
quence star of spectral type F7-8V. 

However, there is currently only transit photometry available 
for WASP-3 from relatively small aperture telescopes (eg the 
IAC 80cm and Keele 60 cm). In this paper we analyse the RISE 
lightcurves to provide improved model fits of the planetary pa- 
rameters and also to demonstrate the potential timing accuracy 
obtainable from the ground using such an instrument. We are 
in the process of targeting several other planetary systems with 
RISE in order to detect TTV signals. 

In §2 we describe the RISE instrument and the observa- 
tion and data reduction procedures. We describe the fitting tech- 
nique used to extract the parameters and uncertainties from the 
lightcurves in §3, and finally in §4 and §5 we present, discuss 
and summarise our results. 
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2. Observations and data reduction 

2.1. The RISE instrument 

RISE was mounted on the LT and commissioned in February 
2008, and was designed to obtain high time precision lightcurves 
of transiting extrasolar planets. Due to transiting systems often 
being relatively bright targets (m v <~ 13), overheads due to 
readout time adversely affect the timing precision we can obtain. 
Therefore, RISE consists of a thermoelectrically cooled frame 
transfer CCD (using an e2v CCD47-20 back illuminated chip) 
which allows exposure times as short as 1 second without dead 
time. For very bright targets (eg HD189733, V ~< 8) this al- 
lows much more data (>10000 images) to be taken during a typ- 
ical 2-3 hour transit than from a traditional CCD with significant 
readout time. 

The CCD has 1024 X 1024 13.0 ftm pixels in the light sen- 
sitive region, giving a field of view of 9.4 x 9.4 arcminutes 
squared. This relatively large field should provide many bright 
sources to be used as reference stars, therefore improving the 
potential photometric accuracy, and allows accurate lightcurves 
to be taken for relatively uncrowded fields previously unattain- 
able with RATCAM - the other imager on the LT with a field of 
view of 4.6 x 4.6 arcminutes squared. 

Furthermore, the RISE instrument uses a single broadband 
filter composed of Schott KG5 and OG515, which results in a 
band pass of 500-700nm. This allows the instrument to collect 
as many photons as possible during an exposure, which again 
should improve the photometric accuracy we can obtain and 
therefore the timing accuracy. For more information about the 
design and performance of RISE, and the LT, see ISteele etaO 
(120081) . 

2.2. RISE photometry 

A full transit of WASP-3 was observed using RISE on the night 
of 18 May 2008, and a partial transit was observed on the night 
of 4 September 2008, correspondin g to epochs E = 250 a nd 309, 
respectively, from the ephemeris of Pol lacco et alj d2008l) . 

T c = 2454143.8503 + 1.846834 £ [HJD]. 

The data from 18 May consist of 4800 x 3s exposures beginning 
at 23:15 UT spanning a 4 hour period starting ~0.8 hours be- 
fore ingress and ending ~0.4 hours after egress, resulting in an 
airmass range of ~ 1.99- 1.02. The 4 September data consist of 
3900 x 3s exposures beginning at 22:30 UT, spanning 3.25 hours 
starting approximately 0.7 hours before ingress and ending just 
before egress - corresponding to an airmass range of ~ 1 .08-2. 10. 

The detector was set to 2x2 binning mode for both transits 
giving a scale of 1 . 1 arcsec/pixel, and the LT was defocussed so 
that the chip was not saturated and well below the non-linearity 
regime resulting in a typical FWHM of ~3-4 pixels (3.3-4.4 arc- 
sec). This also helps to minimise errors associated with flat field- 
ing. Both nights were photometric. 

The images were first de-biased and flat fielded with com- 
bined twilight flats using standard IRAF3 routines. Aperture 
photometry was then performed on the target star and 3 nearby 
companion stars using Pyra£| and the DAOPHOT package us- 

2 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 

3 Pyraf is a product of the Space Telescope Science Institute, which 
is operated by AURA for NASA. 



1.06 - 




0.96 1 1 1 ' ' 1 1 

-0.060 -0.040 -0.020 0.000 0.020 0.040 0.060 

phase 

Fig. 1. RISE lightcurves of WASP-3 taken on 18 May 2008 (up- 
per, offset by 0.05) and 4 Sep 2008 (lower) shown in 1 minute 
bins with the best fit model from the MCMC analysis overplot- 
ted. Residuals are shown offset along with their error bars after 
rescaling so^? . = 1 and to account for red noise. 



ing a 7 pixel aperture (the same companion stars were not used 
in both cases due to different field orientations). The flux of 
WASP-3 was then divided by the sum of the companion stars 
(all checked to be non-variable) to obtain each light curve. Initial 
estimates of the photometric errors were calculated using the 
aperture electron flux, sky and read-noise. The lightcurves were 
then normalised by dividing through with a linear function of 
time fitted to the out-of-transit data, setting the unocculted flux 
of WASP-3 equal to 1, revealing a transit of depth ~1.17%. 
The lightcurves along with their best fit models and residuals 
are shown in Fig. [1] after grouping the data into 1 minute bins. 
The rms of the residuals from the best fit model (see §3.1) are 
~ 4.3 mmag for each lightcurve per 3s exposure. 

3. Results and analysis 

3.1. Determination of system parameters 

In order to determine the system parameters from the photome- 
try, a parameterised model was constructed. This used Kepler's 
Laws and assumed a circular orbit to calculate the normalised 
separation (z) of the planet and star centres as a function of time 
from the stellar mass and radius (M+ and /?*), the planetary mass 
and radius (M p and R p ), the orbital period and inclination (P and 
z), and finally a central transit time for each l ightcurve (7o,„). 
The analytic models of Man del & Agoll d2002l) were then used 
to calculate the stellar flux occulted by the planet from the nor- 
malised separation and the planet/star radius ratio (p) assuming 
the quadratic limb darkening function 

^ = l-fl(l-^)-Ml-^) 2 , 
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where / is the intensity, /j is the cosine of the angle between the 
line of sight and the normal to the stellar surface, and a and b 
are the linear and quadratic limb darkening coefficients, respec- 
tively. 

T he limb darke ning parameters were obtained from the mod- 
els of lClaretl d2000t) . We linearly interpolated the ATLAS tables 
for T eff = 6400 K, log g = 4 25, [Fe/H] = 0.0 and v t = 1.0 km 
s _1 (from Poll acco et al.ll2008l) to obtain limb darkening param- 
eters in both the V and R bands. The average of the parameters 
from each of the bands was then taken to be our theoretical limb 
darkening parameters. 

A Markov Chain Monte Carlo (MCMC) algorithm was then 
used to obtain the best fit parameters a n d their uncertainties 
(see for example, iTegmark et alj l2004t iHolman et al.1 [2006; 



ICollier Cameron et all 2007 ; Winn et al. 20081) . This consists of 
calculating the x 2 fitting statistic, 



XT 



(fj,obs fj,calc) 



where fj, bs is the flux observed at time j, crj is the correspond- 
ing uncertainty and f^aic is the flux calculated from the model 
for time j and for the set of physical parameters described above. 
The best fit parameters and uncertainties are then found by start- 
ing a long chain that explores the probability densities in param- 
eter space by accepting or rejecting subsequent parameter sets 
based on the difference in^ 2 . The new parameter set is accepted 
ifx 2 is lower than that of the previous set, or accepted with prob- 
ability 

-Ay 2 
P = exp(^) 

if x 1 is higher than that of the previous parameter set. New 
parameter sets are chosen by adding a small random amount 
to each of the previously accepted parameters determined by a 
jump function. These jump functions consist of a random value 
selected from a Gaussian distribution with mean and standard 
deviation of 1, scaled by a factor associated with each parameter 
(which ideally should initially be set near to each parameter's 
lcr uncertainty to speed chain convergence). These scale fac- 
tors were then resized from their initial values so that ~25% of 
all new parameter sets were accepted. The reader is referred to 
the aforementioned literature for more details of MCMC fitting 
techniques. 

To obtain reliable estimates of parameters and their uncer- 
tainties it is important that the photometric errors are calculated 
accurately. The photometric errors cr ; are first rescaled so that 
the best fitting model for each lightcurve has a reduced ;f 2 (x 2 ed ) 
of 1 , which required the initial error estimates to be multiplied 

by the factors ^jx 2 rei j i = 1-73 and ^jx 2 re d 2 = 2.04. This indicates 
that the photometric errors are considerably higher than what 
would be expected from photon noise alone. 

It is also vi tal to account for any correlated ( 'red') noise 
in the data (seelPont etalJl2006l: iGillon et al.ll2006l) . Following 
IWinn et al.ld2008l we evaluated the presence of red noise in each 
lightcurve by calculating a factor f3 and again rescaling the <x ; - in 
the x 2 calculation by this value. A value for /3 is determined by 
analysing the residuals from the best fit model of the lightcurves. 
Calcualting the standard deviation of the residuals cr\ and also 
the standard deviation after binning the residuals into M bins of 
N points cr N , one would expect 



cr N 



M 



in the absence of red noise. However, this is u sually larger by 
a factor p. As pointed out in I Winn et al.l d2008l) . J3 only depends 
weakly on the averaging time, therefore we used the average of 
/3 values calculated from a range of averaging times from 10-30 
min to rescale the photometric errors. The values of /3 used were 
1.43 and 1.35 for the full transit and partial transit, respectively. 

To account for any normalisation errors in the lightcurves, 
a further 2 free parameters were added for each transit. These 
were the out-of-transit flux (f oot ,„) and a time gradient (tcmd.n), 
which are vital for TTV measurements as these effect the sym- 
metry of the lightcurve and therefore the central transit times. 
An airmass correction k was tried for the full transit instead of 
the time gradient by multiplying each point by a factor exp(-kz) 
where z is the airmass. This was found to cause an insignificant 
change in the output parameters and so the time gradient was 
used in our final analysis. The chains would not converge when 
an airmass correction was used instead of a time gradient for the 
partial transit. 

A Gaussian prior was added to the x 2 fitting statistic so that 



X 



N 

-z 



(fj,obs ~ fj.cak) 2 (M 



*,J 



■ Mo) 2 



yffi V M - 1 



where M+j is the mass of the star used to calculate each model 
at time j, and Mq and <x M are the stellar mass a nd correspond- 
ing uncertainty as given in iPollacco et al.l ([2008). This is so the 
error in the stellar mass is taken into account when determining 
parameters from the MCMC probability distributions. The stel- 
lar radius cannot be allowed to vary independently of the stellar 
mass as the chains will not converge. Therefore, the stellar radius 
was updated for each choice of M* using the scaling relation 

1 /3 

R* oc MJ~ , whilst P and M p were held fixed at their previously 
determined values, as their uncertainties do not have any signif- 
icant effect on the output probability distributions. 

An initial MCMC analysis was used to estimate the initial 
parameters and jump functions for p, i, 7o„, f 00 t,n an d tcmd,n- 
Five separate chains with 100000 points were then computed 
with the initial free parameters set by adding a 5<x gaussian ran- 
dom to their previously determined best fit values. The first 20% 
of each chain was eliminated to keep the initial conditions from 
influencing the results, and the remaining part of the chains were 
merged to obtain the best fit parameters and uncertainties. These 
were not assumed to be Gaussian and were obtained by setting 
the best fit value to the modal value of the probability distri- 
bution, and the lcr lower (upper) limit as the value where the 
integral of the distribution from the minimum (maximum) value 
to it was equal to 0. 159. To test that the chains had all converged 
to the same parameter space, the Gelman & Rubin statistic was 
then calculated for ea ch of the free parameters (this time assum- 
ing Gaussian errors) (Gelman & Rubin1 [l992h . It was found to 
be less than 0.5% from unity for all parameters, a good sign of 
mixing and convergence. The transit duration {Tj) and impact 
parameter (b) were also calcu lated for each parameter set in th e 
chain using the equations of Seager & Mallen-Ornelasl ([2003), 
and were treated the same as the fitted parameters to determine 
best fit values and uncertainties. This was to account for any cor- 
related errors. 

To check for any errors that may have resulted from a poor 
choice of limb darkening parameters, the above procedure was 
repeated, this time allowing the linear limb darkening parameter 
(a) to vary freely whilst holding the quadratic limb darkening 
parameter (b) fixed at th e theoretical value obtained from the 
models oflClaretl (|2000). This is the same procedure used by 
IWinn et al.l (120081) . and did not cause any significant changes to 



N. P. Gibson et al.: Updated parameters for the transiting exoplanet WASP-3b using RISE 



o 
6 



10 



-10 



-250 



-200 



-150 



-100 
Epoch 



-50 



50 



Fig. 2. Timing resid uals of the RISE tran sits of WASP-3 and 
the epoch given in iPollacco et alJ (I2008I) using the updated 
ephemeris. 



Table 2. Central transit times and uncertainties for the RISE pho- 
tometry. 

Epoch Transit Time [HJD] Uncertaintly (days) 





59 



2454605.55956 
2454714.52210 



0.00035 
0.00036 



tion of Epoch E and period P to the tran sit times from this pap er 
(shown in table |2]i and the ephemeris of Poll acco et al.l (120081) . 

T C (E) = T c (0) + EP, 

where T C (Q) was set to the central transit time from 18 May 
2008 RISE data. The results were T C (Q) = 2454605.55915 ± 
0.00023 [HJD] and P = 1 .846835 ± 0.000002 days. Fig|2]shows 
a plot of the timing residuals of WASP-3 transits using the up- 
dated ephemeris. 



the best fit values, however, it did increase the errors in i and p 
slightly and hence these results were adopted as our final val- 
ues. A further check on the choice of limb darkening parameters 
involved repeating this process by replacing the quadratic limb 
darkening coefficient (b) determined for the combined V+R fil- 
ter by those obtained for the individual V and R filters. This 
caused no significant variation from the previously determined 
results in either case, and thus we conclude that the choice of 
limb darkening parameters had no detrimental effects on our re- 
sults. This was not surprising, as the linear limb darkening pa- 
rameter was much more sensitive to the choice of filter than the 
quadratic parameter. 

4. Results 

4.1. MCMC results 

The parameters and their uncertainties calculated from the 
MCMC fit are shown in Table [1] The planet radius was found 
to be 1.29+o ( i2 Rj, and the planet density 0.82^^ pj - consis- 
tent with that of the discovery paper. Therefore the derived plan- 
etary parameters are still i n agreement with m odels of highly 
irradiated Hot Jupiter s (e.g. lFortney et al.f2007l) as concluded in 
IPollacco et ail d2008l) . 

The inclination is also consistent with that calculated previ- 
ously but it is now known much more accurately with a value 
of 85.06+°-}!? deg - a factor of ~5 increase in precision. However, 
the transit duration was found to be inconsistent with that of the 
discovery paper (~ 4.1<x). Tj is very sensitive to the inclination 
of the system, and we believe the discrepancy arises due to the 
discovery data not being of sufficient quality to constrain it ac- 
curately, and therefore has resulted in an underestimated error in 
the transit duration. If the updated inclination from the RIS E data 
is use d along with the remaining parameters from lPollacco et ail 
(2008), then the calculated is found to be consistent with our 
value. 

4.2. Transit ephemeris 

The RISE transits do not differ significantly from the discovery 
ephemeris however this can be updated using the longer base- 
line. This was found by minimizing x 1 by fitting a linear func- 



5. Summary and Discussion 

We have presented some of the first results from RISE of a full 
and partial transit of the exoplanet system WASP-3. An MCMC 
analysis was used to update t he parameters which were found 
to be consistent with those of IPollacco et al.1 d2008l) . except for 
the calculated value for the transit duration, which is explained 
in terms of the large uncertainty in the previously determined 
inclination. Using the central transit times an updated ephemeris 
was calculated to aid in further observations of this system. 

Although the timing accuracy was not of the order we in- 
tended for RISE (<~ 10s), other transits observed with RISE 
(for TTV exper iments) are produ cing much more impressive 
light curves (see lSteele et al.ll2008l) . This is mainly because the 
photometric errors for each WASP-3 light curve had to be multi- 
plied by the factors ^Jx%d - 1 -73 and 2.04 in our analysis so that 

the reduced x 1 = 1 ■ It is as yet unclear why this is the case, but 
is perhaps due to randomised systematics caused by a combina- 
tion of the moon, changing flat-field structure and poor autogu- 
iding. Further RISE observations currently being made at the LT 
should not suffer from these problems, and we hope to report on 
further transit timing datasets in future papers. 
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Parameter 


Symbol 


Value 


Units 


Planet/Star radius ratio 


P 


1014+ 00010 
u.iuit_o_. 0008 




Orbital inclination 


i 


85.06!]$ 
0-448 + °;!J» 


deg 


Impact parameter 


b 




Transit duration 


T d 


9 7 c^+0.020 
i.. 1 JJ_o.013 


hours 


Transit epoch 


To 


2454605.55915 ± 0.00023 


HJD 


Period 


P 


1.846835 ±0.000002 


days 


Planet radius 


R P 


i 9Q+0.05 


Rj 


Planet mass 


M„ 


1 76 +008 
0-8252 


Mj 


Planet density 


Pp 


Pj 


Planetary surface gravity 


logg P 


o 49+0.O6 

-0.04 
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